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We present a class of models in which dark matter (DM) is a fermionic singlet under the Standard 
Model (SM) gauge group but is charged under a symmetry of flavour that acts as well on the 
SM fermions. Interactions between DM and SM particles are mediated by the scalar fields that 
spontaneously break the flavour symmetry, the so-called flavons. In the case of gauged flavour 
symmetries, the interactions are also mediated by the flavour gauge bosons. We first discuss the 
construction and the generic features of this class of models. Then a concrete example with an 
abelian flavour symmetry is considered. We compute the complementary constraints from the relic 
abundance, direct detection experiments and flavour observables, showing that wide portions of 
the parameter space are still viable. Other possibilities like non-abelian flavour symmetries can be 
analysed within the same framework. 


I. INTRODUCTION 


Establishing the nature of dark matter (DM) is one 
of the fundamental open problems in particle physics 
and cosmology. A weakly interacting massive particle 
(WIMP) is an excellent candidate since (for GeV to TeV 
scale masses) it naturally provides a relic abundance con¬ 
sistent with observations pQ. Since WIMPs must be 
weakly interacting, a likely possibility is that DM is a 
singlet under the Standard Model (SM) gauge group. In 
this case the interactions with the SM particles are trans¬ 
mitted by mediators, i.e. by the “dark sector”. Many 
possibilities for mediators haven been discussed in the lit¬ 
erature, among them Z' models ma, Higgs portal mod¬ 
els mu], z portal [13], pseudoscalar mediation hme 
dark color OMES and models with flavoured DM fT9b28j . 

Flavour models aim at explaining the Yukawa struc¬ 
ture by introducing a flavour symmetry under which the 
SM fermions transform non-trivially. After this symme¬ 
try is spontaneously broken by scalar fields (i.e. flavons), 
the observed fermions masses and mixing angles are gen¬ 
erated as non-renormalisable operators. Froggatt and 
Nielsen proposed the first flavour model by introducing 
an abelian U{ 1) flavour symmetry and heavy vector-like 
quarks [29] . Later on, many other possibilities have been 
studied, among them SU( 2) [30] El], SU( 3) [32[ as well 
as discrete flavour symmetries like A 4 [33] . 

In this article we propose a new mechanism generat¬ 
ing the interactions of DM with the visible sector: DM 
and the SM particles are charged under a flavour sym¬ 
metry and interact with each other only via flavour in¬ 
teractions. 1 This is an appealing possibility since it is 
minimal in the sense that no ad hoc quantum numbers 
must be introduced and DM interactions are described by 


the same dynamics generating fermion masses in the SM. 
Furthermore, relating Dark Matter and flavour models 
provides a handle on the otherwise unspecified flavour¬ 
breaking scale, motivating the possibility of low-energy 
realisations with interesting phenomenology that would 
allow insights into the flavour sector [34]. 


II. GENERAL SETUP 


In this section, we examine the generic features of 
“flavour portal” models that do not depend on the spe¬ 
cific implementation, e.g. the transformation properties 
of SM and DM fields under the flavour symmetry. 

In models a la Froggatt-Nielsen 29l [35] l38l . a new sym¬ 
metry of flavour (Gf), under which the SM fermions are 
charged, is introduced, such that the Yukawa couplings 
of the SM are forbidden at the renormalisable level (with 
the possible exception of the top Yukawa) and fermion 
masses arise once the symmetry is broken as higher di¬ 
mensional operators: 


(mf)ij = a{■ 


\M ) y/2' 


( 1 ) 


Here v = 246 GeV is the vev of the SM Higgs h and 
represents one or more scalar fields, the flavons, whose 
vev breaks the flavour symmetry. M is a cutoff that can 
be interpreted as the mass of vector-like fermions, that 
constitute the UV completion of the model. 2 The expo¬ 
nents n{j are dictated by the transformation properties 

of the fields under Qp , while a{- are unknown coefficients 
originating from the fundamental couplings of SM fields, 
flavons and messengers in the UV-complete theory. They 


1 Note that our framework shares this starting assumption with 
some flavoured DM models, e.g. [2T| [23] 25 |28| . Nevertheless, 
the crucial difference is that we propose either flavour-breaking 
scalars or flavour gauge bosons as the only mediators between 
the dark and the SM sectors. 


2 In Ref. [39] it has been discussed that heavy scalars with the 
quantum numbers of the SM Higgs field can also be used for a 
similar mechanism. Model building details and phenomenology 
of the messenger sector have been discussed in Ref. El ED- 
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are commonly assumed to be 0(1), so that the the ob¬ 
served hierarchies of the Yukawas are solely accounted 
for by the flavour symmetry breaking. This can be done 
by a suitable choice of the flavour charges or representa¬ 
tions of the SM fermions (one can assume that the Higgs 
field is neutral under Gf): provided that e ={(/))/M is a 
small expansion parameter, typically of the order of the 
Cabibbo angle or smaller. 



(a) ( b ) 


A. Gf as a global symmetry 

If the flavour symmetry is global, there are no extra 
gauge bosons associated to it, so that for a DM particle 
which carries flavour charge, the only interactions with 
the visible sector are through flavon fields <fi. In particu¬ 
lar, we are going to consider two chiral fermions xl and 
Xr , which are SM singlets but transform non-trivially un¬ 
der Gf • In perfect analogy with the SM fermions, the DM 
particle acquires a Dirac mass term though the flavour 
symmetry breaking: 



FIG. 1: Feynman diagrams contributing to DM annihilation 
in presence of light flavons (a) and ( b ) and light flavour gauge 
bosons (c). 


that: 



where again b x is an 0(1) coefficient and n x depends on 
the flavour representations/charges of xl and xr- Note 
that in this setup the stability of the DM particles is 
automatically guarantied (at least for what concerns the 
quark sector) by an accidental symmetry related to their 
singlet nature. In fact, due to Lorentz and SU( 3) invari¬ 
ance, DM could only decay by violating baryon number. 
However, like in the SM, baryon number is an accidental 
symmetry of our model, since the dark sector, i.e. the DM 
fields and the DM messengers which constitute the UV 
completion of Eq. ©> is composed only by SM singlets. 
This guarantees the stability of DM at least at the per¬ 
turbative level, since the dark sector fields can not mix 
with quarks and Froggatt-Nielsen messengers (i.e. heavy 
vector-like quarks) because of colour and charge conser¬ 
vation. 

From Eqs. we can obtain the couplings of both 

DM and SM fermions to a dynamical flavon, without 
further specifying the details of the model: 


mj, = k{<f>) . 


( 4 ) 


Since we are not going to specify the details of the scalar 
potential and the symmetry breaking, we take k as a free 
parameter. 

From the above Lagrangian, we can see that our 
setup shares some similarities to the Higgs-portal sce¬ 
nario naHUDHE], in particular the fact that the medi¬ 
ator prefers to couple to heavier fermions, but also the 
possible correlations between the DM annihilation and 
scattering with nuclei. 

In Fig. [l] we show the Feynman diagrams (a) and 
( b ) contributing to DM annihilation in presence of light 
flavons. Notice that the contribution ( b ) requires m x > 
77i0. From the figure, we can immediately infer the para¬ 
metric dependence of the annihilation processes we are 
interested in. Concerning the s-channel processes (a), we 
see that the thermally-averaged annihilation cross section 
mediated by flavons scales as: 


(<4v) 




(4 m\ - ml) 2 + Tim 2 


( 5 ) 


C D n fr ^rf L f R (j) + (n x + 1 ) T ^rXLXR(t> 

w w 

= ^ffLfR<t> + A x XlXR • (3) 

Here, we suppressed flavour indices but one has to keep in 
mind that the couplings to fermions are in general flavour 
changing. In fact, Xf is not diagonal in the same basis 
as the fermion mass matrix m/, as a consequence of the 
flavour dependence of the exponents 77/. The above ex¬ 
pressions can be easily generalised to the case of multiple 
flavons. 

For what concerns the flavon mass, we assume that 
it also arises from the flavour symmetry breaking, such 


where the flavon decay width T 0 depends on Aj and A^. 
From Eq. (|3|, we see that the flavon-mediated process 
preferably involves the heaviest fermions that are kine¬ 
matically accessible and it can be doubly suppressed in 
the case of heavy flavons: both by the propagator and 
by the couplings (that scale as l/ra^). This is the rea¬ 
son why, as we will see in the explicit example of the 
next sections, the correct relic abundance might require 
a resonant enhancement of the annihilation cross section, 
i.e. m x ~ 7770/2. Note also that the annihilation cross 
section in Eq. © depends linearly of the temperature, 
which is a consequence of the velocity-suppression of the 
process. The contribution of the diagram ( b ) of Fig. [l] 




scales like 
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(a» ~ (6) 

in the limit m x and features a p-wave suppression 

as well. 

Elastic DM-nuclei scattering can only be mediated by 
</>, hence being controlled by the same couplings depicted 
in the diagram (a) of Fig. [I] Therefore, the parametric 
dependence of the spin-independent cross section is: 


a 


si 

<t> 


rsj 


\ 2 \ 2 

A X A <f>N 2 
m (t) 


( 7 ) 


Here fi x x is the DM-nucleon reduced mass, and X^x oc 
Xf is the scalar couplings to nucleons uni. 


B. Qf as a local symmetry 


We also consider the possibility that the flavour sym¬ 
metry is gauged, leading to extra gauge boson(s) with 
coupling gp. The phenomenology of the model dras¬ 
tically changes. The reason is that the couplings of 
fermions and dark matter to Z' are completely deter¬ 
mined by the flavour charges (Qx in our abelian exam¬ 
ple) of the particles: 

£ d gFXi^(Q X L p L + Q X r p r)x Z'x 

9F7'f(Qf L PL + Q fR PR)fZ' ll . (8) 


As we can see, the couplings are not suppressed by the 
flavour-breaking scale, as in the case of flavon mediation. 
Furthermore, unlike flavons, flavour gauge bosons do not 
preferably couple to heavy flavours and the DM annihila¬ 
tion, depicted in Fig. [T] (c), is potentially efficient even for 
flavour-breaking scales well above the TeV. In fact, the 
flavour gauge bosons can be substantially lighter than 
flavons if the coupling is weak. In fact, in the abelian 
example above we have mz f = a/2 Qf (</>}• Of course, too 
small values of gp would suppress the annihilation even 
if Z' is light. From the diagram (c) of Fig. [lj we see that: 


(cr Z 'v) 


Qf 


(' m 2 z , — 4m 2 ) 2 + r 2 z ,m 2 z , 




(9) 


where the Z' width Yz> is proportional to g 2 F . Also, as 
we can see the s-channel annihilation is not p-wave sup¬ 
pressed (unlike the flavon-mediation case). 

As in the flavon case the Z’ interactions are in gen¬ 
eral flavour-violating, since quarks of different families 
couple in general differently to the Z ’. This feature is 
going to induce severe constraints on the Z’ mass from 
flavour observables. The direct detection interaction is 
also determined by the gauge boson exchange: 


_SI 

G z , 


Qf^Z'N ,.2 

~~Za V X N‘ 


( 10 ) 


The vector coupling to the nucleons is Xz'N oc gp. 


III. CONSTRUCTION OF AN EXPLICIT 
MODEL 


Let us now consider an explicit example of our general 
idea that DM is communicating with the SM fields via 
flavour interactions. For simplicity we choose a Froggatt- 
Nielsen U(1)f model. 3 We assign flavour charges to the 
SM quarks, Q qi , Q Ui and Q^. The SM scalar doublet 
is neutral under U(1)f and the flavour symmetry is bro¬ 
ken by a single flavon with Q^ = — 1. In addition we 
introduce new fermion which are singlets under the SM 
gauge group but carry U(l)p charges: the dark matter 
particles. 

The SM Yukawas read: 

v?j = <A c<,i+c “ j > 4 = 4 eQqi+Qdj . (u) 

where are assumed to be 0(1) numbers and e = 

( (j))/M . We take e = 0.2 to reproduce the observed 
quarks hierarchies and define the rotations to the fermion 
mass basis as 


(V?)VVE » F, (V?)Vvg = y d , (12) 

where are diagonal matrices. 

Given the hierarchical structure of the Yukawas, the 
rotations are approximated by ratios of the Yukawa en¬ 
tries: 


( vt’% 


(q 


(u,Gf)x 

R 


( u,d ) 

va 

( u,d ) 

4 

( u,d ) 

Vji 

(u,d) 

y\j 


( u,d ) 

%• 

(u,d) ( 

a )j 

n (u,d) 

a j^ { 

n ( u 4) 

a jj 




c Q{u i ,d i )-Q(u F dA) 


(13) 


(14) 


where i > j. From Eq. it follows that the couplings 
to a dynamical flavon defined in Eq. © are given by 


A 


(u,cf) _ (u,cf) 


(Q gi + Q( Uj ,dj)) 


, Qqi+Q(uj,dj ) V 


(15) 


Using the above expressions for the rotations, we can 
then easily estimate the couplings in terms of 

fermion masses and mixing angles. For instance, we have: 


A?2 « (Qqi + Qu 2 )m c {VZ) 12 /(4>), (16) 

ASi «(Q® + O, 1 )m c (VS) 12 /<0). (17) 

As we can see, these couplings are thus expressed in terms 
of physical observables with the residual uncertainty from 
the unknown 0(1) coefficients encoded in the Vp and Vr 


3 In the case the symmetry is global, we should rather consider a 
discrete subgroup Zp C U( 1), in order to avoid massless Nambu- 
Goldstone bosons. For large enough TV, this does not substan¬ 
tially modify the effective theory, such that we can still work in 
terms of charges of a continuous U( 1) [36: . 
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rotations . 4 The DM mass arises as in Eq. Q with the 
exponent given by: 

n x = Qxl ~b Qxr ~ 1* (18) 


Therefore, the model with the global U(1)f is com¬ 
pletely defined - up to 0 ( 1 ) coefficients - by the param¬ 
eters: 


777 , 0 , m x , k = 777,0/(0), (19) 

where a given value of m x can be obtained by suitable 
choices of Q Xl , Q Xr and b x . 

There are few possibilities how U(1)f can be chosen 
such that the measured fermion masses and mixing can 
be reproduced. Out of the possibilities outlined in m, 
here we adopt the following example: 

(Qqii 5 2g 3 ) = (3, 2, 0), 

(Quu 2 u 2 ? Qu 3 ) = (3, 2, 0), 

(S dl , Qd 2 , Qd 3 ) = (4, 2, 2). (20) 


Note that, since for us Q qs = Q Us = 0 , we have A 33 = 0, 
i.e. no unsuppressed coupling of the flavon to the top. As 
a consequence the largest coupling of the flavon is flavour 
violating: 0 t L R c R L . 

In the case U(1)f is local, we additionally have a 
Z' whose couplings are shown in Eq. Given the 

charge assignment above, couplings to light generations 
are larger. Furthermore, once the rotations in Eq. ( 12 ) 
are applied to go to the fermion mass basis, flavour- 
violating couplings arise from Eq. (| 8 |, of the form ^ 
(Qq. — Q qi ) x This induces Z'-mediated FCNC 

at tree-level, setting strong limits on the gauge coupling 
gp for a given Z' mass. 


The above adopted charge assignment is anomalous 
ED, as it is a common feature of U(1)f models that 
successfully reproduce the fermion masses and mixing 
|32 38, 42j|43]. 5 As the usually invoked Green-Schwarz 
mechanism would require a flavour-breaking scale close to 
the Planck scale, we have to assume that anomalies are 
canceled by the unspecified field content of the hidden 
and/or the messenger sector. In particular, if the flavour 
messengers are heavy quarks in vector-like representa¬ 
tions of the SM gauge group, they do not need necessary 
to be vector-like under the flavour symmetry too, leaving 
large freedom to cope with the U(1)f anomalies. 


4 Given the hierarchical structure of the Yukawa, the order of mag¬ 
nitude of the entries of do not change after rotating to the 

mass basis, only the unknown 0(l)s are modified. 

5 However, cf. an anomaly-free solution presented in m and the 
recent study |44M in the context of SU( 5) GUT. 


IV. PHENOMENOLOGICAL ANALYSIS 
A. Flavour Constraints 

As we have seen, both in the global and in the local 
version of our U(1)f example FCNC appear at tree-level. 
Using the bounds on FCNC operators reported in m 
145] . we can estimate the limits on the flavon /Z f mass. 

The strongest limit to the flavon mass comes from the 
(s L d R )(s R d L ) operator, contributing to K — K mixing. 
Integrating out the flavon the above operator is induced 
with the coefficient A^A^i/ttt^. The bounds result: 

A M k : m 0 > Vk x 580 GeV, ( 21 ) 

£k ■ m^ZVkx 0arg(Af 2 A^) x 2.3 TeV, (22) 

where we neglected an overall coefficient, product of fun¬ 
damental 0(1) Yukawa-like couplings, which can still 
conspire to relax the bound to some extent. Although 
the limit from CP violation in K — K mixing is rather 
strong, we see that a mild suppression of the overall 
phase, arg(Af 2 A 2 i) ~ 0.1, is enough to reduce it at the 
level of the bound from A M R - 

In the Z' case, the strongest bounds also come 
from K — K mixing. The leading operator is 
( SL r Y fl dL)(s R 'y fl d R ), whose Wilson coefficient reads: 

^-AQ qiq2 (V^ 12 AQ dl dM)i2, (23) 
777 Z' 

where A Qf 1 f 2 = Qf ± — Q/ 2 . The bounds on the Z' mass 
then result: 


A M k : 

mz> fZ 

( 9F 1 

Vio-b 

| x 210 GeV, 

(24) 

fl< : 

m z > ZZ 

( 9F ) 
\io-b 

| x /arg((V/) 12 (V^) 12 ) 




x 3.3 TeV, 

(25) 


where again we omitted an 0(1) uncertainty due to the 
coefficients entering the rotations in Eqs. (fl3{ [l4|). The 
stringent bound from e R can be relaxed at the level of 
the CP conserving limit if arg ((U/)i 2 (Ur)i 2 ) ~ 0 . 01 . 

The fields that constitute the UV completion of 
Froggatt-Nielsen models (vector-like quarks or heavy 
scalars) can give contributions to the FCNC that are 
larger than those mediated by flavon exchanges |34], even 
though they enter at the one-loop level. The reason is 
that the flavon couplings are proportional to the fermion 
masses, suppressing processes involving light generations. 
Adopting the model-independent approach of [34] . we 
find that, in our model, the strongest bound in the 
hypothesis of suppressed phases to the messenger scale 
comes from D — D mixing and it can be translated to a 
limit on the flavon mass: 

^40 ^ Cd x k x 2.3 TeV, 


(26) 
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where Cp parameterises a product of unknown 0(1) co¬ 
efficients. In presence of 0(1) CP-violating phases, we 
obtain the following bound from 

m (f) >CKxkx27 TeV. (27) 

As it will be clear from the discussion in the following 
subsection, the phenomenologically interesting region of 
the parameter space would be excluded by this limit, 
unless we assume that the overall phase in Ck can reduce 
it by about one order of magnitude. 


B. Relic abundance and Direct Detection 


We have implemented our model in the MicrOmegas 
code m in order to obtain an accurate numerical cal¬ 
culation of the relic abundance and the direct detection 
cross section. In the case of a global U(1)f, we work 
with the free parameters m^, ra x , k = /{(/>) and set 

the 0(1) coefficients to unity. In the gauged case, we 
have in addition mz' and gp and the dependence on un¬ 
known coefficients (in the fermion rotation matrices) is 
much milder. We require the relic abundance not to over¬ 
shoot the Planck measurement [I], taking as a conserva¬ 
tive limit Udm^ 2 < 0.13. 

Global U(1)f case. The results for the global case are 
shown in the first plot of Fig.[2j As we can see, besides the 
flavon resonance, there is a wide region with m x > 
where efficient annihilation is provided by diagram ( b ) 
of Fig. [l] and FCNC constraints are satisfied. We display 
here only the bound of Eq. (21) from flavon exchange, as¬ 


suming a mild suppression of the CP-violating phases at 
the level discussed in the previous subsection. We neglect 
the more stringent bounds of Eq. (26) as they rely on 


the additional assumption of a weakly-coupled messenger 
sector at the scale M. We just notice that they would 
exclude the region where yy —>> provides the correct 


relic density, leaving the resonance at m^ 


ty /2 


the only viable solution, unless a mild suppression comes 
from the coefficient Cd- This would be also the effect 


of the bound from e^, cf. Eq. (22), in presence of 0(1) 
phases. The points shaded in grey at low values of 
do not fulfil the LUX constraints m , shown in the sec¬ 
ond plot. As we can see, most of the parameter space 
has good prospects to be tested by next generation di¬ 
rect searches experiments and the only points that might 
be hidden under the neutrino background correspond to 
resonant annihilating DM with m x > 1 TeV. 

Local U(l)p case. In the local case, we find that the 
only viable possibility to fulfil the relic density bounds 
relies on resonant Z' exchange, m x « mz'/ 2, as a con¬ 


sequence of the stringent FCNC constraint of Eq. (25). 
Hence, we varied gp and mz f and scanned m x around the 
resonant condition. The result is shown in the right plot 
of Fig. [2j We see that for a given ra^, too small values 
of gp cannot give the correct relic density, even on the 
resonance, while there is an upper bound on gp from the 
FCNC bound of Eq. (25) and the limit given by direct 


detection experiments is comparably much weaker. As a 
result, only a narrow region, depicted in green, is still vi¬ 
able. The best way to probe this surviving region seems 
to rely on an increased sensitivity of FCNC constraints, 
especially in the D — D and K — K systems. 


C. Collider signals 

There are in principle very interesting signals that can 
be searched for at the LHC. In fact, the largest branching 
ratio of the flavon is </>—)► £c, tc, which is around 60%, but 
also the bb or bs channels are sizeable (while the BR to XX 
is very suppressed because it is at threshold or below due 
to relic density constraints), leading for example to final 
states like pp tctc or pp —» bbbb. However, the produc¬ 
tion cross section of 0 is small, at the level of 0.1(10 -3 ) fb 
for LHC@14 TeV for flavons with = 500 (1000) GeV, 
values already at the edge of the flavour constraints. A 
similar situation happens in the Z' case, where flavour 
bounds make any LHC signal very challenging. 

Finally, let us mention that Higgs-flavon mixing can 
be induced by quartic terms in the scalar potential such 
as as well as by loops of SM fermions. We 

checked that in the former case the mixing scales like 
(u/(</>)) 2 , while it is suppressed by a further factor v/{(/>) 
in the latter case, due to an additional flavon-fermion 
vertex. As a consequence, any effect on the properties of 
the observed Higgs will be suppressed by at least a factor 
(u/(</>)) 4 so that the mixing will affect the h decay widths 
at most at the level of few percent for </> satisfying the 
flavour bounds discussed above. We leave a more detailed 
discussion of this interesting aspect of low-energy flavour 
models to future research. 


V. CONCLUSIONS AND OUTLOOK 

In this article we discussed the possibility that the DM 
particle is a singlet under the SM gauge group but car¬ 
ries flavour charges, such that it can communicate with 
the SM only through flavour gauge bosons or flavour¬ 
breaking scalars. We discussed the general features of 
such flavour portal scenario. We studied an explicit ex¬ 
ample with an abelian 17(1) flavour symmetry. For this 
model the relevant phenomenological constraints are relic 
density, direct detection and flavour observables. These 
measurements often constrain the model to be close to 
the resonant regime m x ~ m^^z' /2, but still future direct 
detection experiments will be able to probe an important 
part of the parameter space of the model. 

Extensions of these kinds of models to the leptonic 
sector are straightforward. Without entering the details, 
here we just notice that, for the sake of dark matter sta¬ 
bility, mixing of the x singlets with neutrinos must be 
avoided, e.g. by enforcing a conserved symmetry (a nat¬ 
ural choice would be the lepton number itself) or by a 
suitable construction of the UV completion of the model. 
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FIG. 2: Left : points with Qdm h 2 < 0.13 in the (ra x , ra^) plane for different values of k. The dashed lines represents the 
corresponding lower bounds on m<$> from FCNC constraints. Centre : nucleon-DM scattering cross section scaled by the actual 
DM density £ • crsi (with £ = Ll x h 2 / 0.11) for the same points as before. Right : (#f, M Z ') plane in the local U(1)f case; only 
the green band is allowed by all data. 


While we illustrated our concept for a specific model, 
the basic features of flavour portals models as discussed 
in Sec. [TT] are generic. On the other hand, the details of 
the phenomenology will depend of the specific realisation. 
Therefore it will be interesting to extend the above dis¬ 
cussion to different charges assignments or other flavour 
groups, in particular non-abelian symmetries. 
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